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ABSTRACT 

We have performed Smoothed Particle Magnetohydrodynamics (SPMHD) simulations 
demonstrating the production of collimated jets during collapse of IMq molecular cloud 
cores to form the 'first hydrostatic core' in low mass star formation. Recently a number of 
candidate first core objects have been observed, including L1448 IRS2E, L1451-mm and Per 
Bolo 58, although it is not yet clear that these are first hydrostatic cores. Recent observations 
of Per Bolo 58 in particular appear to show collimated, bipolar outflows which are incon- 
sistent with previous theoretical expectations. We show that low mass first cores can indeed 
produce tightly collimated jets (opening angles < 10°) with speeds of ~ 2-7 km/s, consistent 
with some of the observed candidates. We have also demonstrated, for the first time, that such 
phenomena can be successfully captured in SPMHD simulations. 
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1 INTRODUCTION 

Majestic, collimated jets are a defining hallmark of the the star 
formation process, observed during the earliest stages of star for- 
mation when the nascent young protostar remains deeply em- 
bedded in its parental molecular cloud. Magnetic fields are the 
thought to provide the main mechanism for launching such out- 
flows during protost ellar c ollapse — via either a centrifugal 'fling' 
Blandford & Pavnj [l982h or a magnetic pressure-driven 'spring' 
Lvnden-B elL i2003h — though the details are only just beginning 



to be und erstood via numeri cal models of the star formation pro- 
cess (e.g. ISeifried et al .^ 2011). The earliest models of protostel- 
lar collapse bylLarsonl (1969), performed in spherical symmetry, 
showed that gravitational collapse evolves through two distinct 
'stair phases: First, the gas becomes optically thick to radiation 
— trapped due to dust opacity — leading to adiabatic heating and 
the formation of the "first hydrostatic core". This ph ase persists 
for a relatively short ti n ie (10^-10^ years. e.^.JBoss & Yorke..l995.: 



Masunaga et al.lll998l : iMasunaga & Inutsuk^ l2000l : iTomida et aP 
20ld : lBatell201lh before the dissociation of molecular hydrogen 



leads to the onset of a secondary collapse phase to form the sec- 
ond, or stellar, core. Recent numerical models have followed the 
first and second collapse in three dimensions with incre asing physi - 
cal realism: beginning with barotr opic hvdrodvnamics (|B atel 1 998h . 
magnetohydrodynamics (MHD) (Baneriee & Pudritz 2006), resis- 
tive magnetohydrodynamics (Machida et al. 2006, 2008) an d ra- 
diative transfer in the flux-limited diffusion approximation (JBatd 
12011.) . 

Magnetohydrodynamical calculations of either the first or both 
stages of collapse of isolated molecular cloud cores tend to show 



relatively low velocity (v ^ 2 kn i/s), wide-angled outflows pro- 
duced during the first core phase (Tomisaka 2002; Machida et al.| 
2006, 2008; Hennebelle & Froma ng 2008; Commerc on et al. 2010 ; 
iBurzle et al.ll201lh with faster, high velocity (v > 30 km/s) and 
well collimated outflow s only produced during the second collapse 
to form the protostar (JMachida et all 120061; iBaneriee & Pudrit a 
l2006l : lMachida et al.l l2008Y The production of collimated outflows 
from the second collapse phase is in line with observed proto- 
stellar jets which are very well collima ted and contain veloci- 
ties of up to several hundred km/s (e.g. iReipurth & BallvluOOll ; 
iBallv. Reipurth & DavisI 120071) . consistent with the escape veloc- 
ity very close to the protostellar surface. Slower and less colli- 
mated outflows produced in simulations are generall y associated 
with the observed wide angled molecular outflows (IRicher et al.l 
l200Ql : lArce eTaDl2007l : lHatchen et al.ll2007l) . 

Recent observations have given tantalising hints of objects de- 
tected while still in the earliest, first hydrostatic core phase: LI 448 
IRS2E (IChen et al.l20ld). L1451-mm dPineda et al.l2 01 1). and Per- 
Bolo 58 (lEnoch et al.|[201Ql : iDunham et all 1201 ih . While L1451- 
mm shows hints of a slow (~ 1.3-1.7 km/s for L1451-mm, see 
IPineda et al.l201 l'). wide outflow consistent with simulations of the 
first core phase, [Dunham et al. (201 1) report a slow (v > 2.9 km/s) 
but well-collimated outflow (opening angles of ^ 8°) associated 
with Per-Bolo 58. This is in contrast to simulations which tend to 
show poorly collimated outflows during the first core phase. 

In this Letter, we report on simulations, performed in ideal 
magnetohydrodynamics, of the collapse of an idealised IMq 
molecular cloud core to the first core phase. These demonstrate that 
highly collimated jets with velocities of up to ^ 7 km/s and open- 
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ing angles of < 10° may indeed be produced from the first core, 
provided the degree of ionisation is sufficiently high. 



2 METHODS 

We solve the equations of self-gravitating, ideal MHD given by 
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where p is the density, v is the velocity, P is the hydrodynamic 
pressure, B is the magnetic field, is the gravitational potential 
and po is the permeability of free space. We solve these equa- 
tions using a standard Smoothed Particle Magnetohydrodynamics 
(SPMHD) scheme, evolving B/p as t he magnetic field variable 
(Eq. [3]), using the _B0rve, Omang & Trulsen (2001.) source-term 
approach for stability, and with artificial viscosity and resistivity 
terms added to capture shocks and magnetic discontinuities, respec- 
tively JPrice & Monaghan 2005). We use time dependent artificial 
viscosity and resistivity parameters as described in Price (i2012h . 
here using aAv G [0.1, 1] and as G [0,0.1]. In particular, we 
do not employ the Euler p otentials approach, as used in p revious 
star formation simulations ("Price & Bate"2007', '2008', '2009), which 
means that there is no restriction on the geometry or winding of the 
field in our simulations. Instead, we control magnetic divergence 
errors by employing a v ersion of the h yperbolic divergence clean- 
ing scheme proposed bv lPrice & Monag han (2005i [),but revise d and 
substantially improved as described in Tricco & Price (^2012^. We 
find that this approach reduces the divergence errors in the collaps- 
ing core by at least two orders of magnitude, substantially reducing 
the non-conservation errors in momentum that previously resulted 
in our inability to evolve beyond the collapse phase with a B /p- 
based approach (see e.g. Pr ice & Federratlv ^2010': Price 2011). It 
also improves on the approach employed by Biirzle et al. (2011), 
where artificial resistivity alone was used to control the divergence 
errors. We found artificial resistivity alone insufficient for stable 
and accurate long-term evolution of the jet/outflows found here. 
In this paper we use a simple barotropic equation of state 
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where Cs is the isothermal sound speed, pc 
pd — 10~^°g/cm^. Particles initially in the external medium are 
assigned a higher Cs corresponding to the higher temperature. Sink 
particles (IBate et al.lll995h of radius 5 AU are inserted once the 
peak density exceeds 10~^°g/cm^, meaning that in this Letter we 
restrict our study to the first core phase only. 

The initial conditions are a IM© dense, cold spherical, uni- 
form density and slowly rotating core in pressure equilibrium 
with a warm, low density ambient medium. The core has radius 
i?c = 4 X lO^^cm (2.7 x 10^ AU), giving an initial density of 
po — 7.4 X 10~^^ g/cm^ and a gravitational free-fall time of 
iff = 2.4 X 10^ yrs. We use Cs = 2.2 x 10^ cm s"\ The 
core is placed inside a larger, cubic domain with a side length of 
8 X 10^^ cm and a 30 : 1 density ratio between the core and the 



ambient medium, in pressure equilibrium, giving a sound speed 
in the external medium of Cs, medium = 1.2 x 10^ cm s~^, or 
Cs, medium — 4:.2GM/R SO that the self- gravity of the external 
medium is irrelevant. For simplicity we use periodic but non- self- 
gravitating boundary conditions on the global domain. The core is 
set in solid body rotation with Q = 1.77 x 10~^^rad s~^, corre- 
sponding to a ratio of rotational to gravitational energy /3r ^ 0.005 
and r^tff = 0.14. We use 1 x 10^ equal mass SPH particles in the 
core, with the density ratio giving 4.8 x 10^ particles in the sur- 
rounding medium. We also performed calculations using 3 x 10^ 
particles in the core. R esolving the Jeans length according to the 
iBate & Burkertl(ll997h criterion would require ^ 3 x 10^ particles. 
The magnetic field is initially uniform in the z-direction, with 
strength Bo characterised by the parameter p, specifying the mass- 
to-magnetic flux ratio (M/^) in units of the critical value for a 
unifo rm spherical cloud (e.g. iMestell 1 19991 : iMac Low & KlessenI 
l2004l) . 
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where ci is a p a ramete r determined numerically by 
iMouschovias & Spitzej (Il976h to be ci :^ 0.53. We have 
performed simulations over a range of magnetic field strengths 
(p = 20, 10, 7.5, 5, 4 and 3), but the main calculations we show 
employ p = 5, corresponding to Bq = 163/iG in physical units 
and an initial plasma /3 (ratio of gas to magnetic pressure) of 3.3 
in the core. It should be noted that the initially imposed vertical 
field is extremely weak compared to the 10-300 mG fields that are 
wound up to produce the jet. 

The setup described abov e is otherwise identical to that em- 
ployed bv lPrice & Bate! (l2007h (hereafter .PB 07.) . except that a bug 
later discovered in the code meant that, although IpBOTI stated that 
their cores were in pressure equilibrium, in fact the calculations 
in that paper were performed with zero pressure in the external 
medium. The main effect of this is that the collapse time is slightly 
longer and that the gas pressure at a given time is slightly lower in 
the collapsing core. We have here performed calculations both with 
and without an external confining pressure, finding similarly well- 
collimated jets, but at a slightly lower initial field strength when the 
external pressure is absent, consistent with the plasma [3 being the 
main parameter controlling the launch/collimation of the outflows. 
Jets and outflows were not produced in the IpBOTI calculations be- 
cause of the restrictions on the field geometry imposed by the Eu- 
ler potentials formulation, meaning that the winding of the toroidal 
field that launches the jet (c.f. Sec. [3]) was not captured. 



3 RESULTS 

Fig. [T] shows snapshots of column density from the p = 5 calcula- 
tion taken every 0.02 free-fall times (490 years) from t/tf[ = 1.03- 
1.13, shortly after the formation of the first core. The wind up of 
toroidal magnetic fields in the inner disc leads to the launch of a 
strong, well-collimated (< 10° opening angle) jet, with outflow ve- 
locities of ^ 5 km/s. The jet itself can be seen to entrain a slower, 
wider outflow and at later times (rightmost panels) shows distinct 
'kinks' or 'wiggles' as a result of material entrained in the helical 
magnetic field expanding in the z direction. At ^ 3, 000 years after 
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Figure 1. Column density rendering showing the launch and propagation of the first-core jet and the entrained outflow from our simulations of a collapsing 
IMq core, showing snapshots every 0.02 free-fall times (490 years ). The well-collimated jet at the centre of the flow expands outwards at 3-7 km/s with an 
opening angle of < 10°, similar to the outflow recently observed bv lDunham et al.l ( 12011]) from the candidate first hydrostatic core Per Bolo 58. 



first core formation (rightmost panel) the outflow has already ex- 
panded to several thousand AU, beyond our initial core radius and 
well into the surrounding medium. The jet continues to be driven 
essentially until all of the collapsing material has been used up (this 
occurs at t/t^^ ?^ 1.2, after which material from the initially hot ex- 
ternal medium starts to be accreted, though the outflow continues 
to be driven as long as mass is supplied). Although a flattened, disc- 
like object appears surrounding the central object, orbital velocities 
in the midplane are sub-Keplerian by a factor of ^^3-4. 

Fig.|2]shows the accreted mass (i.e., the mass of the sink parti- 
cle, red short-dashed line) and ejected mass (defined as all material 
with a spherical Vr > 0.1 km/s, green long-dashed line) as a func- 



tion of time in the simulation, as well as the sum of these (solid 
line). The outflow is remarkably efficient — with up to 40% of the 
initial material in the core ejected. The lower panel shows the max- 
imum velocity of jet material (solid line) together with the mean 
velocity of all particles with Vr > 0.1 km/s, indicating maximum 
speeds of 5-7 km/s in the jet with a mean velocity of outflowing 
material around 2 km/s. 

Fig.[3]shows a rendering of the magnetic field in the simulation 
at t/tff = 1.1 (showing the field on each SPMHD particle drawn 
with an opacity proportional to the field strength). The field near 
the protostar is strong (^-^ 100 mG) and tightly wound in a toroidal 
geometry, which is responsible for the high degree of collimation. 
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Figure 2. Cumulative total accreted and ejected mass as a function of time 
(top panel, as indicated), together with the maximum velocity of ejected 
material (bottom panel, solid/black line) and the mean velocity of all par- 
ticles with Vr > 0.1 km/s (bottom panel, red/dashed line). The outflow is 
remarkably efficient, expelling 40% of accreted material by the end of the 
simulation (halted once the complete initial core mass has been accreted) at 
maximum speeds of 5-7 km/s (bottom panel), with a mean velocity of all 
outflowing material ~ 2 km/s. 



Animations of the field evolution show clearly the growth of this 
"magnetic tower" flow (Lvnden-Bell 2003). 

We find similar jets are produced for a range of magnetic field 
strengths, specifically for 4 < /x < 10 in our setup (Fig. |4]), cor- 
responding to initial field strengths of Bq ^ 80-200/iG and initial 
plasma /3's in the range 2-13. The /i = 10 calculation produces a 
transient jet but then only a wide angled, poorly collimated wind, 
while simulations with /x > 4 tend to induce strong magnetic brak- 
ing of the central object. Weaker field (/x < 20) calculations are 
complicated by the fact that secondary fragmentation tends to oc- 
cur in the present setup when a barotropic equation of state is used 
(c.f. lMrzleetal. 2011) . 



4 DISCUSSION 

In this Letter we have demonstrated the production of highly col- 
limated jets with velocities of up to 7 km/s during the first hy- 
drostatic core phase of star formation. The features are similar to 
the collimated jets fou nd in other sirnulations, either during the 
second collapse ph ase (Tom isakall2002l : iBaneriee & Pu dritz 200d: 
iMachida et al.ll2008 ), at much higher masses dSeifried et al...201lh 
or at low mass but with lower velocities dCiardi & Hennebellel 
|201Q|) . 

Firstly, our simulations serve to demonstrate, for the first time, 
the production of stable, well collimated jets with a Smoothed 
Particl e Magnetohydrodynamics (SPMHD) technique. Biirzle et al.l 
(l201lh demonstrat ed that a standard artifi cial-resistivity based 
SPMHD approach jPrice & MonaghanluOOS ) could be used to sim- 
ulate protostellar outflows, without the restrictions on magnetic 
field geometr y associated with earlier Euler-potentials-based ap- 
proaches fe.g. IPrice & Batell200l 12008.) . We find that using arti- 
ficial resistivity alone to control the divergence error is insufficient 
for long term, stable evolution of the jets found here and that a high 
resistivity is required simply to control the divergence error. Our 
implementation of a new divergence cleaning technique (described 





Figure 3. Rendering of the magnetic field geometry (left) and magnetic cur- 
rent (right) in the /x = 5 simulation, showing the strongly wound, toroidal 
magnetic field and strong currents that launch the jet. Wiggles and kinks 
seen in the column density (Fig.[TJ can also be seen in the field geometry, 
indicating that this material is entrained by the magnetic field. Magnetic 
field strengths in the jet are ~ 10-300 mG. 



in detail in lTricco & Pricell2012h means that we have been able to 
follow the stable growth of the jet to several thousand AU — until 
the entire molecular cloud core has been accreted — without requir- 
ing a high resistivity, which we found could suppress the formation 
of the collimated flow. 

The jets we find are similar to those produced in other cal- 
culations, both protostellar and in accretion discs more generally 
(see review by IPudritz et al .1120 07), driven by a 'magnetic tower' 
of tightly wound toroidal magnetic field expanding at a veloc- 
ity related to the orbital velocity at the footpoint of the tower 
(iLvnden-B elll2003h . Wiggles similar to those seen here (e.g. Fig.[T]) 
and a corkscrew-like morphology are a common feature arising 
from 3D jet simulations (Ouyed et al. j2003l) . 

We find the main requirements for the production of colli- 
mated outflows during the first core to be weak, but not too weak, 
initial magnetic fields (specifically, 4 < /x < 10 for the setup 
described here), rota tional motion, and relatively low resistivity. 
iMachida et al.l (120061) also note that the jets produced in their (sec- 
ond core) simulations are sensitive to resistivity, noting that jets 
always occur in their calculations employin g ideal MHD, but not 
always in their resistive MHD calculations. ICiardi & Hennebellel 
(2010) find similar, collimated outflows from the first core in their 
ideal MHD simulations, albeit with lower outflow speeds (v ~ 1.2- 
1.8 km/s) and dependent on the angle between the magnetic field 
and the rotation axis. 

Considering the high physical resistivity thought to be present 
in molecular clouds at these densiti es (e.g. Nakano et al. 2002), the 
use of ideal MHD, both here and in lCiardi & Hennebell^ (l2010l) . is 
extremely wfzrealistic, and one would not exp ect jets to be produced 
from the first core on this basis. However, the lDunham et al.l (1201 ih 
observations do suggest a collimated, bipolar outflow from a first 
core candidate. One possibility is that first-core jets are still pos- 
sible at high resistivity (77 > 10^° cm^/s) but require higher mag- 
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Figure 4. Effect of varying the initial magnetic field strength. With weaker 
fields (/i = 10, left panel) an initial fast jet is launched but detaches, leaving 
a transitory 'knot' in a wide angled, poorly coUimated wind. At higher field 
strengths the jet is more tightly coUimated (/i = 5 and /i = 4, middle and 
right panels, respectively). 



netic field strengths in the initial conditions or only occur at higher 
density. Preliminary investigations in this regard suggest that this 
may indeed be the case. The other possibility is that the conduc- 
tivity in Per-B olo 58 is indeed higher than the estimates given by 
iNakano et al.l (2002), perhaps due to thermal ionisation following 
first core formation. Further simulations should be able to distin- 
guish between these two possibilities. Measurements of the mag- 
netic field strength and/or geometry in Per-Bolo 58 would also pro- 
vide a critical constraint on simulation models. 

Our simulations suggest that objects such as Per-Bolo 58, 
showing a coUimated, bipolar outflow with a characteristic veloc- 
ity of ^ 3 km/s, may indeed be viable first core candidates. The 
outflow in L1451-mm is les s well coUimated, but remarkably effi- 
cient — IPineda et al.l (120111) note that the amount of mass needed 
to keep material at 560 AU with a velocity of 1.3 km/s is ^ 0.53 
M0, almost twice the mass observed in the dense core itself. This 
is consistent w ith the high efficiency found here (Fig. |2l see also 
ICiardi & Hen nebelle 2010), where up to 40% of the mass in our 
original core is ejected (and more if mass is continually supplied 
from an external reservoir). On the other hand , speeds of ^^ 25 km/s 
reported for L1448 IRS2E (IChen et aLl2010h would seem inconsis- 
tent with an object in the first core phase, though it should be noted 
that velocities > 9 km/s would not be reached in the simulations 
shown here due to the 5 AU size of the central sink particle. 
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